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Abstract The forms of Euler and Gibbs-Duhem relations discussed in thermodynamics 
of extensive systems are shown to hold also for nonextensive systems with long-range 
interactions with a novel interpretation of the entities appearing therein. In this way, the 
principles underlying Tsallis' scaling relations in equilibrium nonextensive thermo- 
statistics are clarified. 
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Nonextensive systems with long-range interactions reside outside of the realm of 
ordinary thermostatistics [1-3]. Their macroscopic quantities do not scale linearly with 
respect to system size and therefore the standard thermodynamic limit cannot naively be 
employed. It is important to notice that they are "small systems" since their elements 
always feel the effects of the boundaries through the long-range interactions. It seems that 
these systems have been regarded for a long time as anomalous ones that are of less 
importance from the pragmatic viewpoint of thermostatistics. The situation has changed, 
however. For example, there are active experimental and theoretical studies on non- 
neutral plasmas [4,5]. Also, development of numerical analysis of self-gravitating 
systems [6] stimulates investigation of traditional gravothermodynamics in a new 
perspective [7]. One may also add to these the nanometric systems of great current 
interest [8-12] in view of smallness of the systems. In particular, Hill and his co-worker 
[8,9] have developed for small systems a modified version of the conventional 
thermodynamics, which preserves the Euler relation for the internal energy but violates 
the Gibbs-Duhem relation. In the present paper, we present a different theory, in which 
the forms of the Euler as well as Gibbs-Duhem relations are retained but with different 
interpretations of the terms appearing therein. 

Let us imagine a generic model with a long-range interaction. As an example, consider 
a spin model defined on a J-dimensional hypercubic lattice whose Hamiltonian reads 
H = -J 52 <7, o j I r", where cr, is the spin on the ith site, 7 is a positive coupling 

constant, r tj is the distance between the ith and jth sites, and a(>0) is a constant 
characterizing the property of the interaction range. In the continuum approximation, the 
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internal energy per spin can be estimated to be [13,14] U I N N* , where N is the total 
number of spins. N* is a factor given as follows: 



N~ 



const (a > d) 

XnN (a = d) (1) 
N l - a,d (0<a<d), 



in the large-A^ limit. This factor highlights the nonextensive feature of the system with 
a < d and, in fact, universally appears in various systems with long-range interactions 
[15,16]. The asymptotic behavior of N* in Eq. (1) is determined by the property of the 
power-law decay of the interaction. 

The scaling discussion in Ref. [13] is inspired by the nonextensive internal energy, 
U oc N N* . There, the Gibbs free energy of a nonextensive system is assumed to be of the 
same form as that of an ordinary extensive system and to satisfy the following scaling 
property: 



G U T S P V 

■ + T7^-> (2) 



NN NN" N N N* N' 

where T, S, P, and V are the temperature, the entropy, the pressure, and the volume, 
respectively. (The magnetic energy may also be added, if necessary.) This equation 
means that the entropy (as well as the volume) scales with and, therefore, is extensive, 
whereas the ordinary intensive quantities, the temperature and the pressure, scale with ./V* 
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[13,14,16]. This turns out to be consistent with the Monte Carlo results reported in Ref. 
[17] (see also Ref. [18]). 

In what follows, first we reconsider the Gibbs free energy and its scaling property with 
the viewpoint of retaining the forms of the Euler and Gibbs-Duhem relations, but with 
altered interpretations of the entities entering therein. Here, we are concerned with the 
concept of homogeneity and, therefore, we limit ourselves to the strongly nonextensive 
case 



Let us start our discussion by considering the internal energy, U(S,V,N), of a 
nonextensive system with a long-range interaction, which scales as 



where a is a function of the entropy and volume per particle, s = S I N and v = V I N, 
respectively. Therefore, we have 



< a < d. 



(3) 



U(S, V, N) = u(s, v)NN*, 



(4) 



U(XS, XV, XN) = X' 



2-a/d 



U(S, V, N). 



(5) 



Taking the derivatives of the both sides with respect to A, we have 



d(XS) d(XV) 



dU r dU 
— — S+— — 



v+ 



du 



N = (2-a/d)X 



l-ald 



u. 



(6) 
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Putting 



1 



l-ald 



1 = 



(7) 



\2-al d ) 



we obtain the Euler relation 



U = TS-PV + LlN, 



(8) 



where T = dU/dS, P = -dU I dV , and the chemical potential fi = dU I d N. Note that 
in the discussion of the extensive systems, one sets X equal to unity. The choice in Eq. 
(7) is a reflection of the nature of the problem of long-range interactions. Thus, we obtain 
our first result that the Euler relation for nonextensive systems is formally equivalent to 
the ordinary Euler relation [19] for extensive systems. This means that the form of the 
Gibbs-Duhem relation also remains unchanged. From Eqs. (6) and (8), it follows that 



which is, in fact, formally equivalent to the ordinary Gibbs-Duhem relation for extensive 
systems. Eq. (8) yields the thermodynamic potential, G- jlN , and therefore Eq. (2) is 

justified. Since G <^ N N* , the chemical potential scales with /V*, similarly to T and P, 



SdT-VdP + Ndn = 0, 



(9) 
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indicating the necessity of incorporating the fluctuations in the system. Consequently, we 
obtain the following scaling relations: 



These relations should be compared with Eqs. (3.20)-(3.22) in Ref. [19], in which the 
prefactor, 1/ N* , is absent. The new Gibbs-Duhem relation presented in Eq. (9) shows 
that there is a new relationship between the parameters defined by Eqs. (10)-(12). With 
respect to these variables, we may define other thermodynamic quantities such as specific 
heat, compressibility and so on with similar definitions of magnetization and 
susceptibility. This justifies the manipulation employed in Ref. [17], in which Monte 
Carlo calculations of long-range interacting Ising-like systems are performed. 

Finally, let us illustrate the above development by taking S(U, V) to be of the form 




(10) 




(11) 




(12) 



S = cU^ V 



(13) 



with constant c, £, , and 77 . Then, extensivity of S leads to the following condition on £ 
and 77: 
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(2-a/<0§ + 7/ = l. 



(14) 



Then, we obtain the following expressions for the thermodynamic parameters in Eqs. 
(10)-(12): 



T 1 U ( S\ 



N* %NN\Nj 



(15) 



N* % NN* 



N 



-i 



P n U fV , 

- ' : (16) 



4 = ^^- (17) 
/Y £ 7Y7Y 



In these forms, the scaling properties of the quantities and variables are manifestly 
appreciated. 

In conclusion, we have established the Euler and Gibbs-Duhem relations for 
equilibrium nonextensive thermostatistics of the systems with long-range interactions. 
We have seen how this approach can clarify the principles underlying Tsallis' scaling 
relations. It is of crucial importance to notice that in the present discussion we are 
concerned with equilibrium. It is generally considered that the large system-size limit and 
the long-time limit do not commute with each other in nonextensive systems. If the long- 
time limit is taken first, one observes nonequilibrium stationary states, which survive for 
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very long periods [20]. Currently, there is an ongoing investigation of understanding 
physical properties of such states based on a generalization of conventional Boltzmann- 
Gibbs statistical mechanics, termed nonextensive statistical mechanics [21-23], which 
employs Tsallis' nonadditive entropy [24]. In this respect, attention should be paid to the 
recent discussion [25] that Tsallis' scaling relations are valid only in equilibrium states, 
and not in nonequilibrium stationary states. 

Note Added. The Euler relation in nonextensive systems has been considered in Ref. [26], 
without reference to the scaling properties. Also, the authors of Ref. [27] discuss the 
Gibbs-Duhem relation in terms of Tsallis' nonadditive entropy. Therefore, there are no 
overlaps between the present work and these ones. 
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